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Bloch  Wave  Approach  to  the  Optics  of  Crystals: 
The  Role  of  Spatial  Dispersion 
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Dipartimento  di  Fisica  del  Politecnico  di  Torino 
Corso  Duca  degli  Abruzzi  24,  10129  Torino,  Italy 
Fax:  + 39-011-564  7399;  email:  ponti@polito.it 

Abstract 

In  this  work  we  give  new  results  in  the  framework  of  a research  on  the  macroscopic  optical 
properties  of  crystals  whose  period  p is  shorter  than  the  light  wavelength  but  large  on  the 
mQlecular  scale,  so  that  they  can  be  treated  in  the  framework  of  a continuum  theory.  More 
precisely,  we  assume  that  the  optical  properties  of  the  crystal  are  fully  defined  by  a perfectly 
periodic  function  e(r ),  representing  the  local  permittivity  of  the  medium  (assumption  of  locality). 
To  define  the  effective  permittivity  tensor  e we  use  a method  based  on  the  analysis  of  the  normal 
modes  for  the  electromagnetic  field  propagating  within  of  the  periodic  medium  (Bloch  waves). 

1.  Introduction 

A perfect  crystal  is  generally  treated  in  optics  as  a homogeneous  medium,  despite  the  fact  that 
in  crystals  the  time  harmonic  solutions  of  Maxwell  equations  (i.e.  the  normal  modes  of  the 
electromagnetic  field)  are  Bloch  waves,  whereas  in  homogeneous  media  they  are  plane  waves. 
We  recall  that  a Bloch  wave  can  be  considered  as  a plane  wave  whose  amplitude  is  a periodic 
function  of  r,  or  equivalently  as  a superposition  of  plane  waves,  obtained  by  expanding  the 
amplitude  in  a Fourier  series.  The  dielectric  tensor  e of  the  homogeneous  ( macroscopic ) model 
for  crystals  is  defined  by  considering  as  a starting  point  the  Bloch  wave,  and  neglecting  the 
space  dependence  of  its  amplitude.  This  means  that  only  the  plane  wave  component  with  q = 0, 
that  defines  the  macroscopic  field  [1,  2,  3],  is  taken  into  account,  and  considered  as  the  best 
‘plane  wave  approximation’  of  the  actual  Bloch  wave.  The  homogeneous  medium,  where  this 
plane  wave  can  freely  propagate,  is  the  effective  macroscopic  medium.  Its  tensor  e is  therefore 
implicitly  defined  by  the  equation 


Do  ~ £ (£o^o)  ’ (1) 

where  Eg  define  the  zeroth  order  Fourier  components  of  the  Bloch  wave  and  eo  is  the  vacuum 
permittivity. 

Here  we  discuss  the  application  of  the  Bloch  wave  method  to  the  study  of  a short  period 
crystal  whose  optical  properties  are  assumed  as  perfectly  known  and  fully  defined  by  a peri- 
odic tensor  field  e(r).  This  field  represents  a mesoscopic  model  for  the  actual  crystal,  where 
its  molecular  ( microscopic ) structure  is  neglected,  and  the  assumption  of  locality  is  explicitly 
made.  Despite  this  last  assumption,  the  spatial  dispersion  plays  a main  role  in  our  theory, 
since  e explicitly  depends  on  the  wave  vector  k of  the  internal  plane  wave,  i.e.  e — I(k).  The 
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tensor  field  e(r)  is  well  known  for  many  periodic  liquid  crystal  (LC)  phases  and  for  artificially 
made  structures  (as  for  instance  Reusch  piles).  Our  problem  reduces  therefore  to  the  search 
of  equations  relating  e(k)  to  e(r).  Fully  analytic  expressions  for  e(k)  will  be  given,  that  allow 
for  a quantitative  analysis  of  the  approximations  required  to  define  homogeneous  models  for 
non-homogeneous  and  periodic  media.  This  is  the  main  motivation  of  our  research.  The  found 
equations  will  be  applied  to  the  particularly  interesting  and  simple  case  of  cholesteric  liquid 
crystals. 

2.  Theory 

The  effective  dielectric  tensor  e of  the  macroscopic  model  is  expressed  as  a function  of  the  Fourier 
components  of  the  tensor  field  e(r): 

£(r>  = J2  e?  exP  (w*0-  (2) 

9 

The  component  of  order  zero  is  the  space  average  of  e(f),  i.e.  — e,  and  constitutes  a first 
rough  approximation  for  e.  The  other  approximations  are  obtained  by  adding  to  e terms  having 
the  general  structure: 


£Qi  G-qi  £92  (91  +92)  (91  +...+9at  _ ) ) £Qn  1 (3) 

where 

Gg  = ^[{k  + q)2I-(k  + q)(k  + q)-e\  \ (4) 

I is  the  3x3  identity  matrix,  (k  + q)  (k  + q)  is  a dyadic  product,  and  the  vectors  qn  satisfy  the 
following  relations: 


N' 

9i  + - + = 0,  \tN'<N.  (5) 

n=l 

The  terms  given  by  Eq.  (3)  can  be  interpreted  as  the  effect  of  the  multiple  scattering  within  the 
periodic  structure,  and  the  equation  (5)  states  that  only  the  forward  scattering  with  multiplicity 
N >2  gives  a contribution  to  e.  The  formal  expression  of  e can  be  written  as: 

OO 

£~£+  ^2  e9i^— 9l£92^— (9l+92)-"^— (9i  + -.+9N-i)e-(9i+-+9N-i)-  (®) 

N— 2 9i  9W-1 

The  dependence  of  e on  k comes  from  the  tensors  G$,  as  shown  by  Eq.  (4).  The  dominant 
contribution  to  e is  given  in  general  by  the  terms  with  N = 2,  i.e.  by  the  double  scattering 
(two-photon  scattering),  a fact  that  allows  to  greatly  simplify  Eq.  (6).  A further  simplification 
is  obtained  by  expanding  Gg  in  power  series  of  p/X,  where  p is  an  average  period  of  the  crys- 
tal lattice,  since  only  the  first  terms  of  such  expansion  are  expected  to  give  a non-negligible 
contribution  to  e(k).  This  gives  * 


e = E + e^°)  -f 


*2) 


+ ... 


(7) 
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For  the  simple  case  of  crystals  periodic  in  only  one  direction,  say  X3  (ID  crystals),  the  vectors 
q depend  on  a single  index  r,  and: 

qr  = rqx3,  i},  = 2_^eTG'_r 'er; 

r^O 


where: 
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(r  = 0, ±1,  ±2, ...),  q = 27 r/p,  n = k/k 0,  £33  = -^(£13^1  + £23^2))  and  only  the  two-photon 
scattering  has  been  considered.  For  chiral  structures,  the  most  important  term  is  the  one 
containing  and  scaling  as  p/ A,  since  it  is  related  to  the  optical  activity  of  the  structure. 
Interestingly,  it  is  identically  zero  for  light  propagating  along  the  periodicity  axis  of  ID  crystals, 
i.e.  for  n — 713X3,  as  clearly  shown  by  the  expression  of  Gr*\  Eq.  (8).  The  higher  order  terms 
give  in  general  small  corrections  to  the  preceding  ones,  but  in  some  particular  case  they  can 
give  qualitatively  new  effects.  This  actually  occurs  in  cubic  crystals  (and  in  particular  in  the 
blue  phases),  because  the  term  scaling  as  (p/A)2  gives  a small  anisotropy  to  the  structure.  This 
also  occurs  along  the  periodicity  axis  of  ID  crystals,  whose  optical  rotation  is  related  to  a term 
scaling  as  (p/A)3,  since  the  term  scaling  as  p/A  is  zero.  In  fact,  for  real  e(r)  only  the  terms 
scaling  as  (p/A)m  with  odd  m-values  give  a contribution  to  the  optical  activity. 


3.  Discussion 

(i)  The  method  developed  here  allows  to  easily  obtain  e(k)  in  all  the  cases  where  the  function 
e(r)  has  only  few  non-zero  Fourier  components,  a fact  that  actually  occurs  in  many  liquid  crystal 
phases.  The  simplest  and  most  important  one  is  the  choelsteric  phase,  where  only  the  Fourier 
components  eT  with  r = 0,  ±1  are  different  from  zero.  For  such  crystals: 
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and  £^  = 0,  where  e0  = (e0  + £e)/2,  £e  = ee,  e'  = (ee  - £0)2(2£0  - n\  - v%)/ 8 and  e"  = 
i(ee  — £0)2(2e0  — nf  — n2)/(4e0).  The  pedices  0 and  e refer  to  ordinary  and  extraordinary, 
respectively.  We  recall  that  cholesteric  liquid  crystals  are  locally  uniaxial,  with  the  optic  axis 
orthogonal  to  X3,  and  that  the  pitch  of  their  helical  structure  is  equal  to  2p,  where  p is  the 
lattice  period.  The  Eqs.  (9)  show  that  effective  medium  is  still  uniaxial,  but  with  the  optic  axis 
parallel  to  X3.  The  term  scaling  as  (p/A)  is  identically  zero  for  any  direction  of  the  light  beam: 
the  optical  activity  of  cholesteric  is  therefore  given  by  the  term  scaling  as  (p/A)3.  This  unusual 
property  constitutes  a further  demonstration  of  the  unique  optical  properties  of  the  cholesteric 
phase. 

One  may  notice  that  the  equations  given  here  can  straightforwardly  give  the  expression  of 
e up  to  terms  scaling  as  (p/A)3,  whereas  the  standard  effective- medium  theories  generally  only 
give  terms  scaling  as  (p/A)°  and  (p/A)1. 

(it)  The  presence  of  ^-dependent  terms  give  non-easy  problems  for  the  search  of  plane  wave 
solutions  and  for  the  boundary  conditions.  In  fact  for  any  given  direction  of  k the  dispersion 
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relation  becomes  a polynomial  that  could  have,  in  principle,  any  number  of  roots.  However,  only 
four  solutions,  ±&i  and  ±&2,  corresponding  to  plane  waves  with  different  polarization  states, 
have  physical  meaning.  The  difficulty  for  the  boundary  conditions  is  due  to  the  fact  that  the 
terms  of  e(k)  depending  on  the  m-power  of  k come  from  the  space  derivatives  of  order  m of  the 
electric  field.  The  constitutive  equations,  written  in  the  Landau  form,  are  in  fact: 

% = ^if  + lijkd/dxk  + Xijkid*  / dxkdxi  + ...;  p = p0. 

Interestingly,  the  Bloch  wave  method  gives  the  actual  expression  of  e(k)  without  any  use  of 
higher  order  tensors.  In  the  presence  of  derivative-depending  terms,  the  usual  conditions  of 
continuity  for  the  tangential  components  of  E and  H axe  no  more  valid.  The  problem  has  been 
at  least  partially  solved  only  for  the  simple  case  of  first-order  derivatives  [5],  by  writing  the 
constitutive  equations  in  the  Post  or  Tellegen  form  [6]. 

(in)  The  limits  of  validity  of  the  macroscopic  model  have  been  carefully  tested.  For  achiral 
crystals,  the  approximation  l(k)  ~ e(0),  where  all  the  fc-dependent  terms  are  omitted,  is  in 
most  cases  good  enough  up  to  p/X  « 0.1.  The  full  expression  of  e(k)  allows  to  extend  these 
limits  up  to  p/X  « 0.5,  except  for  the  optical  properties  directly  related  to  the  spatial  dispersion 
(as,  for  instance,  the  rotatory  power).  Such  properties  require  a more  careful  analysis,  that  has 
been  partially  given  on  the  basis  of  a different  approach  in  Ref.  [7].  The  limit  p/X  m 0.5  is  due 
to  the  fact  that  for  higher  p- values  the  periodic  medium  can  give  Bragg  diffraction,  where  the 
dispersion  curves  show  forbidden  bands.  Obviously,  homogeneous  models  are  not  able  to  account 
for  the  Bragg  diffraction  bands,  even  though  they  could  still  give  approximated  expressions  for 
some  optical  properties  at  higher  p- values. 
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